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Abstract 

We review major progress on the modeling of electric dipole emission from rapidly 
spinning tiny dust grains, including polycyclic aromatic hydrocarbons (PAHs). We 
begin by summarizing the original model of spinning dust proposed by Draine and 
Lazarian and recent theoretical results improving the Draine and Lazarian model. The 
review is focused on important physical effects that were disregarded in earlier studies 
for the sake of simplicity and recently accounted for by us, including grain wobbling 
due to internal relaxation, impulsive excitation by single-ion collisions, the triaxiality 
of grain shape, charge fluctuations, and the turbulent nature of astrophysical environ- 
ments. Implications of the spinning dust for constraining physical properties of ultra- 
small dust grains and environment conditions are discussed. We discuss the alignment 
of tiny dust grains and possibility of polarized spinning dust emission. Suggestions for 
constraining the alignment of tiny grains and polarization of spinning dust are also 
discussed. 

1 Introduction 

Diffuse Galactic microwave emission carries important information on the fundamental prop- 
erties of the interstellar medium, but it also interferes with Cosmic Microwave Background 
(CMB) experiments (see Bouchet et al. [4]; Tegmark et al. [71]). Precision cosmology with 
Wilkinson Microwave Anisotropy Probe ( WMAP) and Planck satellite requires a good model 
of the microwave foreground emission to allow for reliable subtraction of Galactic contami- 
nation from the CMB radiation. 

The discovery of an anomalous microwave emission (hereafter AME) in the range from 
10 — 100 GHz illustrates well the treacherous nature of dust. Until very recently it has 
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been thought that there are three major components of the diffuse Galactic foreground: 
synchrotron emission, free-free radiation from plasma (thermal bremsstrahlung) and thermal 
emission from dust. In the microwave range, the latter is subdominant, leaving essentially 
two components. However, it is exactly in this range that an anomalous emission component 
was reported (Kogut et al. [44] and [45]). In the paper by de Oliveira-Costa et al. [11] this 
emission was nicknamed "Foreground X" , which properly reflects its mysterious nature. This 
component is spatially correlated with 100 fim thermal dust emission, but its intensity is 
much higher than one would expect by directly extrapolating the thermal dust emission 
spectrum to the microwave range. 

An early explanation for AME was proposed by Draine & Lazatian [20] and [19] (hereafter 
DL98 model), where it was identified as electric dipole emission from very small grains 
(mostly containing polycyclic aromatic hydrocarbons-PAHs) that spin rapidly due to several 
processes, including gas-grain interactions and dust infrared emission. Although spinning 
dust emission had been discussed previously (see Erickson [27]; Ferrera & Dettmar [29]), 
Draine & Lazarian were the first to include the variety of excitation and damping processes 
that are relevant for very small grains. 

While the DL98 model appears to be in general agreement with observations (see [53]; 
[30]), it did not account for some important effects, namely, the non-sphericity of grain 
shapes, internal relaxation within grain, and transient spin-up due to ion collisions. 

This induced more recent work in order to improve the original DL98 model. The recent 
papers include (Ali-Hai'moud et al. [1]; Hoang et al. [36]; Ysard [79]; Hoang et al. [41]; Silsbee 
et al. [70]). In this paper we review both the original DL98 model and the ways that it has 
been improved recently. We focus on the improvement of dynamics of PAHs and important 
physical effects associated with these ultrasmall grains. Recent reviews of the subject include 
Draine & Lazarian [22] and [57], and Lazarian & Finkbeiner [54]. 

In §2 we briefly present the history of AME and discuss the original DL98 model including 
their basic assumptions. §3 presents our principal results improving the DL98 model from 
Hoang et al. [36] and [41]. From §4 to §6 we review on grain rotational dynamics and discuss 
our general approach to calculate power spectrum of spinning dust emission, grain angular 
momentum distribution and emissivity for PAHs of arbitrary shapes. In §7, we discuss 
the implications of spinning dust for constraining physical parameters of PAHs as well as 
environment conditions. The possibility of polarization of spinning dust and its constraint 
is discussed in §8. A summary of the present review is given in §9. 

2 The original DL98 model 

2.1 Anomalous microwave emission and PAHs 

The emission spectrum of diffuse interstellar dust was mostly obtained by the InfraRed 
Astronomy Satellite (IRAS) and infrared spectrometers on the COsmic Background Explorer 
(COBE) and on the InfraRed Telescope in Space (IRTS). The emission at short wavelength 
(A < 50/im) arises from transiently heated ultrasmall grains (e.g., PAHs). These grains have 
such a small heat capacity that the absorption of a single ultraviolet (UV) starlight photon 
(~ 6 eV) raises their temperature to T vib > 200 K. Typically, these grains have less than 300 
atoms and can be viewed as large molecules rather than dust particles. They are, however, 
sufficiently numerous to account for most of the prominent 2175A absorption feature, and 
for ~ 35% of the total starlight absorption (see e.g. Li & Draine [61]). 
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The thermal (vibrational) emissivity of these grains is thought to be negligible at low 
frequency, because they spend most of their time cold and only emit most of their energy 
when they are hot. These ultrasmall grains (PAHs) are invoked in the DL98 model to account 
for the anomalous microwave emission (AME) that was measured in observations. 

The first detection of anomalous dust-correlated emission by COBE (Kogut et al. [44], 
[45]) was quickly followed by detections in the data sets from Saskatoon (de Oliveira-Costa 
et al. [9]), OVRO (Leitch et al. [60]), the 19 GHz survey (de Oliveira-Costa et al. [13]), and 
Tenerife ([12]). Initially, AME was identified as thermal bremsstrahlung from ionized gas 
correlated with dust (Kogut et al. [44]) and presumably produced by photoionized cloud rims 
(McCullough et al. [65]). This idea was scrutinized in Draine & Lazarian [20] and criticized 
on energetic grounds. Poor correlation of Ha with 100 /im emission also argued against 
the free-free explanation (McCullough et al. [65]). These arguments are summarized in 
[22]. Later [10] used Wisconsin H- Alpha Mapper (WHAM) survey data and established that 
the free-free emission "is about an order of magnitude below Foreground X over the entire 
range of frequencies and latitudes where it is detected". The authors concluded that the 
Foreground X cannot be explained as the free-free emission. Additional evidence supporting 
this conclusion has come from a study at 5, 8 and 10 GHz by Finkbeiner et al. [33] of several 
dark clouds and HII regions, two of which show a significantly rising spectrum from 5 to 10 
GHz. 

The recent Wilkinson Microwave Anisotropy Probe (WMAP) data were used to claim a 
lower limit of 5% for the spinning dust fraction at 23 GHz (Bennet et al. [2]). However, other 
models of spinning dust are not ruled out by the WMAP data, and in fact fit reasonably 
well. Finkbeiner [30] performed a fit to WMAP data using a CMB template, a free-free 
template (based on Hct-correlated emission plus hot gas emission near the Galactic center), 
a soft synchrotron template traced by the 408 MHz map, a thermal dust extrapolation 
(Finkbeiner et al. [31]) and a spinning dust template consisting of dust column density times 
Tj. This fit results in excellent x 2 /dof values (1.6, 1.09, 1.08, 1.05, 1.08) at (23, 33, 41, 61, 
94) GHz and a reasonable spectral shape for the average spinning dust spectrum. 

This WMAP analysis alone does not rule out the Bennet et al. [2] hypothesis of hard 
synchrotron emission, but when combined with the Green Bank Galactic Plane survey data 
(Langston et al. [48]) at 8 and 14 GHz, spinning dust appears to provide a much better fit 
than hard synchrotron (Finkbeiner et al. [32]). 

Spinning dust emission has recently been reported in a wide range of astrophysical en- 
vironments, including general ISM (Gold et al. [34]; [35]; Planck Collaboration [67]), star 
forming regions in the nearby galaxy NGC 6946 (Scaife et al. [69]; [73]), and Persus and 
Ophiuchus clouds (Casassus et al. [7]; Tibbs et al. [72]). Early Planck results have been 
interpreted as showing a microwave emission excess from spinning dust in the Magellanic 
Clouds (Bot et al. [3]; Planck Collaboration [67]). 

2.2 Basic assumptions 

(i) The smallest PAH particles of a few Angstroms are expected to be planar. The grain 
size a is defined as the radius of an equivalent sphere of the same mass. PAHs are assumed 
to be planar, disk-like with height L and radius R for a < and spherical for a > a<i- The 
value a2 = 6 A is adopted. 

(ii) PAHs usually have electric dipole moment n arising from asymmetric polar molecules 
or substructures {intrinsic dipole moment) and from the asymmetric distribution of grain 
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charge. The latter is shown to be less important. 

(iii) The grain spins around its symmetry axis ai with angular momentum J parallel to 
ai and J is isotropically oriented in space. 

(iv) For a fixed angular momentum, the spinning grain emits electric dipole radiation at 
a unique frequency mode u, which is equal to the rotational frequency, i.e., v = oo/2n. 

(v) A grain in the gas experiences collisions with neutral atoms and ions, interacts with 
passing ions (plasma-grain interactions), emits infrared photons following UV absorption, 
and emits electric dipole radiation. All these processes result in the damping and excitation 
of grain rotation, i.e., they change grain angular momentum J and velocity ui. 

(vi) Due to the excitation of various aforementioned processes, the grain angular velocity 
randomly fluctuates and its distribution can be approximated as the Maxwellian distribution 
function / Mw (w). 

(vi) The total emissivity per H atom of the electric dipole radiation from spinning dust 
at the frequency v is given by 



where rin is the density of H nuclei, is the electric dipole moment perpendicular to the 
rotation axis, and dn/da is the grain size distribution function with a in the range from a m ; n 



3 Improved Model of Spinning Dust Emission 

Ali-Haimoud et al. [1] revisited the spinning dust model and presented an analytic solution 
of the Fokker-Planck (FP) equation that describes the rotational excitation of a spherical 
grain if the discrete nature of impulses from single-ion collisions can be neglected. 

Hoang et al. [36] (hereafter HDL10) improved the DL98 model by accounting for a number 
of physical effects. The main modifications in their improved model of spinning dust emission 
are as follows. 

(i) Disk-like grains rotate with their grain symmetry axis a! not perfectly aligned with 
angular momentum J. The disaligned rotation of J with ai causes the wobbling of the grain 
principal axes with respect to J due to internal thermal fluctuations. 

(ii) The power spectrum of a freely spinning grain is obtained using Fourier transform. 

(iii) Distribution function for grain angular momentum J and velocity u are obtained 
exactly using the Langevin equation (LE) for the evolution of J in an inertial coordinate 
system. 

(iv) The limiting cases of fast internal relaxation and no internal relaxation are both 
considered for calculations of the angular momentum distribution and emissivity of spinning 
dust. 

(v) Infrequent collisions of single ions which deposit an angular momentum larger than 
the grain angular momentum prior to the collision are treated as Poisson-distributed events. 

The wobbling disk-like grain has anisotropic rotational damping and excitation. Such an 
anisotropy can increase the peak emissivity by a factor ~ 2, and increases the peak frequency 
by a factor 1.4 — 1.8, compared to the results from the DL98 model. 

The effects of grain wobbling on electric dipole emission was independently studied in 
Silsbee et al. [70] using the FP equation approach, but they disregarded the transient spin-up 
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by infrequent single-ion collisions, and considered two limiting cases of dust grain tempera- 
ture Td — > and Td — > oo. 

Further improvements of the DL98 model were performed in Hoang et al. [41], where a 
couple of additional effects were taken into account: 

(i) emission from very small grains of triaxial ellipsoid (irregular) shape with the principal 
moments of inertia Ji > I2 > I3. 

(ii) effects of the orientation of dipole moment n within grain body for different regimes 
of internal thermal fluctuations. 

(iii) effects of compressible turbulence on the spinning dust emission. 

The work found that a freely rotating irregular grain with a given angular momentum 
radiates at multiple frequency modes. The resulting spinning dust spectrum has peak fre- 
quency and emissivity increasing with the degree of grain shape irregularity, which is defined 
by Ii : I2 : I3. Considering the transient heating of grains by UV photons, the study found 
that the spinning dust emissivity for the case of strong thermal fluctuations is less sensi- 
tive to the orientation of /x than in the case of weak thermal fluctuations. In addition, the 
emission in a turbulent medium increases by a factor from 1.2-1.4 relative to that in a uni- 
form medium, as sonic Mach number M s increases from 2-7. The latter Mach numbers are 
relevant to cold phases of the ISM (see Hoang et al. [41] for more details). 

4 Grain Rotational Configuration and Power Spectrum 

A discussion of the basic physical processes involved in spinning dust can be found in the 
review by Yacine Ali-Haimoud, which can be found in the same volume. There, the use 
of Fokker-Planck equation for describing grain dynamics is discussed. Here we discuss our 
numerical approach based Fourier transform and the Langevin equation, which exhibits a 
number of advantages to the FP equation when numerical studies of grain dynamics are 
performed and arbitrary shape of PAHs is considered. We summarize a general approach to 
find the spinning dust emissivity from grains of triaxial ellipsoid shape with l\ > I2 > I3 
subject to fast internal relaxation. 

4.1 Torque- free motion and internal relaxation 

The dynamics of a triaxial (irregular) grain is more complicated than that of a disk-like grain 
with I 2 = I3. Indeed, in addition to the precession of the axis of major inertia &i around J 
as in the disk-like grain, the axis a x wobbles rapidly, resulting in the variation of the angle 
6 between ai and J (see Figure 1). 

To describe the torque- free motion of an irregular grain having a rotational energy E TOt , 
the conserved quantities are taken, including the angular momentum J, and a dimensionless 
parameter that characterizes the deviation of the grain rotational energy from its minimum 
value, 



Q = 




(2) 



The orientation of the triaxial grain in the lab system is completely described by three 
Euler angles -0, and 9 (see e.g., Hoang et al. [41]). Following [74], we define the total 
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Figure 1: Rotational configuration of a triaxial ellipsoid characterized by three principal axes 
ai, &2 and a3 in the inertial coordinate system xyz. Grain angular momentum J is conserved 
in the absence of external torques and directed along z-axis. The torque-free motion of the 
triaxial grain comprises the rotation around the axis of major inertia ai, the precession of 
ai around J, and the wobbling of ai with respect to J. The dipole moment /x, which is fixed 
to grain body, moves together with the grain and thus radiates electric dipole emission. 



number of states s in phase space for q ranging from 1 to q as 



s = 1 



2 

7T JO 



dip 



h(h - hq) + h(h - h) cos 2 j) 



2„/,l 1 /2 



(3) 



where 



COS 



h{hq-h) 
h(h-h) 



nl/2 



(4) 



for q > q sp and ipi = 7r/2 for q < q sp , with g sp = I\jl<i being the separatrix between the two 
regimes. 

The Intramolecular Vibrational- Rotational Energy Transfer process (IVRET) due to im- 
perfect elasticity occurs on a timescale 10 -2 s, for a grain of a few Angstroms (Purcell [68]), 
which is shorter than the IR emission time. So, when the vibrational energy decreases due to 
IR emission, as long as the Vibrational- Rotational (V-R) energy exchange exists, interactions 
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between vibrational and rotational systems maintain a thermal equilibrium, i.e., T rot « T virj . 
As a result, the LTE distribution function of rotational energy reads (hereafter VRE regime; 
see Lazarian & Roberge [58]): 

/vreO, J) oc exp (- ^"f ) « exp f- ^ ) . (5) 

Substituting i? rot as a function of J and q from Equation (2) into Equation (5), the 
distribution function for the rotational energy becomes 

/ VEE ( S ,J) = ^exp(-^iy, (6) 
where A is a normalization constant such that Jq 1 /vre(s, J)ds = 1. 



4.2 Power spectrum of a freely spinning grain 

Consider a grain with a dipole moment fi fixed in the grain body rotating with an angular 
momentum J. If the grain only spins around its symmetry axis, then the rotating dipole 
moment emits radiation at a unique frequency v equal to the rotational frequency, i.e., v = 
u)/27r (see DL98). The power spectrum for this case is simply a Delta function 8{v — lv/27t) 
with a unique frequency mode. 

For an irregular grain of triaxial ellipsoid shape, the grain rotational dynamics is more 
complicated. In general, one can also obtain analytical expressions for power spectrum, but 
it is rather tedious. To find the power spectrum of a freely rotating irregular grain, Hoang 
et al. [36] and [41]) have employed a more simple, brute force approach based on the Fourier 
transform approach. First, they represent the dipole moment n in an inertial coordinate 
system, and then compute its second derivative. We obtain 

3 

A» = 53/i i a i , (7) 

j=i 

where //j are components of y, along principal second derivative of with 

respect to time, and i — 1, 2 and 3. 

The instantaneous emission power by the rotating dipole moment is equal to 

P cA (J,q,t) = ^ 2 . (8) 

The power spectrum is then obtained from the Fourier transform (FT) for the components 
of /i. For example, the amplitude of jl x at the frequency is defined as 

/+oo 
il x (t)exp(-i27rv k t)dt, (9) 
-oo 

where k denotes the frequency mode. The emission power at the positive frequency v\. is 
given by 

PedM <?) = Afe + fj-lk + Aifc)> ( 10 ) 
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where the factor 2 arises from the positive/negative frequency symmetry of the Fourier 
spectrum. To reduce the spectral leakage in the FT, we convolve the time-dependent function 
jj, with the Blackman- Harris window function (Harris 1978). The power spectrum then needs 
to be corrected for the power loss due to the window function. 

The total emission power from all frequency modes for a given J and q then becomes 

P cd (J,q) = £P cd , fe (J,g) = \ \] dt , (11) 

where T is the integration time. 1 

Figure 2 presents normalized power spectra (squared amplitude of Fourier transforms), 
|FT(/i X)J/ )| 2 /max(|FT( / u x )| 2 ) and |FT( / u 2 )| 2 /max(|FT( / u x )| 2 ) for the components jl x (or jl y ) 
and fl z , for a freely rotating irregular grain having the ratio of moments of inertia I\ : I<i : 
J 3 = 1 : 0.6 : 0.5 and for various q. Circles and triangles indicated with m and n denote peaks 
of the power spectrum for oscillating components of jl x (or jl y ) and jl z , respectively. The 
horizontal axis is the angular frequency of emission modes normalized over the frequency of 
emission when the grain spins around its shortest axis. 

Multiple frequency modes are observed in the power spectra of the irregular grain, but 
in Figure 2 we show only the modes with power no less than 10~ 3 the maximum value. 
One can see that for the case with large q = 1.6, the modes with oj/(J/I\) > 1 have power 
increasing, while the modes with oj/(J/I\) < 1 have power decreasing. It indicates that if 
grain rotational energy is increased so that the grain spends a significant fraction of time 
rotating with large q, then the grain should radiates larger rotational emission. 

Although one should not expect the analytical expression of power spectrum for the 
triaxial grain, the frequency modes can be approximately found. Indeed, for q < I1/I2, we 
found that power spectra for jl x (or jly) have angular frequency modes 

u m ^(4>)+m(\iP\), (12) 

where the bracket denotes the averaging value over time, and m = 0,±1,±2... denote the 
order of the mode. The frequency modes for jl z are given by 

u n = n(\j>\), (13) 

where n is integer and n > 1. 

In the following, the emission modes induced by the oscillation of \i x or fi y , which lie 
in the xy plane, perpendicular to J, are called in-plane modes, and those induced by the 
oscillation of \l z in the direction perpendicular to the xy plane, are called out-of-plane modes. 
The order of mode is denoted by m and n, respectively. Figure 2 also shows that the emission 
power for out-of-plane modes u n is negligible compared to the power emitted by in-plane 
modes u m . 

Emission power spectra are numerically calculated to find Uk and P e d,k, as functions of 
J and q, for the various ratio of moments of inertia I\ : I2 '■ I3. The obtained data will be 
used later to compute spinning dust emissivity. 



: This is the result of Parseval's Theorem. 
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Figure 2: Normalized power spectrum of a torque-free rotating irregular grain with Ii : I 2 : 
I 3 = 1 : 0.6 : 0.5 for the different values of q — 1.05,1.60 (i.e. q < q sp = h/h) and q = 
1.81 > g sp . The components of |FT(/i a; )| 2 /max(|FT(/i a; )| 2 )) (or |FT(/i y )| 2 / max(|FT(/i 2 .)| 2 )) 
are indicated by circles, while the components of |FT(/i^)| 2 / maxdFT^)! 2 ) are indicated 
by triangles. Orders of in-plane modes m and out-of-plane modes n are indicated, and case 
1 (/ii = /i/v^) of n orientation is assumed. Figure reproduced from Hoang et al. [41]. 



5 Grain Angular Momentum Distribution: Langevin 
Equation 

5.1 Langevin Equation 

To find the exact distribution function for grain angular momentum J, Hoang et al. [36] and 
Hoang et al.[41] proposed a numerical approach based on the Langevin equation. Basically, 
they numerically solved the Langevin equation describing the evolution of three components 
of J in an inertial coordinate system. They read 

dJi = Aidt + ^Biidqi, for % = x, y, z, (14) 

where dqi are random Gaussian variables with (dqf) = dt, Ai = (AJi/At) and Bn = 
((AJj) 2 /At) are damping and diffusion coefficients defined in the inertial coordinate sys- 
tem. Detailed expressions of these coefficients can be found in Hoang et al. [36] and Hoang 
et al. [41]. 

For an irregular grain, to simplify calculations, we adopt the Ai and Bn for a disk-like 
grain obtained in HDL10. Following DL98b and HDL10, the disk-like grain has the radius 
R and thickness L = 3.35 A, and the ratio of moments of inertia along and perpendicular 
to the grain symmetry axis h = I\\/I±- Thereby, the effect of nonaxisymmetry on Ai and 
Bn is ignored, and we only examine the effect of grain wobbling resulting from the grain 
triaxiality. 

1 /2 

In dimensionless units, J' = J//||0;t,|| with cjt,|| = (2A;B^ g as/-f||) being the thermal 
angular velocity of the grain along the grain symmetry axis, and t' = £/ r H,|h Equation (14) 
becomes 



dJ[ = A\dH + JB'M 



(15) 
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where (dq 2 ) = dt' and 



V 2 J' 3 

A', = -pi— (16) 

' gas.eff ° ' ed,eff 

B'u = 27^-^11, (17) 



where 



/ T gas,eff _ tot, || 

gas ' cff 7S,|| ~ cos 2 ^ + 7Hsin 2 e' 1 1 

^tot,±TH,|| / Tcd.cff , lrA 

7H = p 7^ eff = , (19) 

-Ttot,||TH,± T H ,|| 

where th,|| and th,_l are rotational damping times due to gas of purely hydrogen atom for 
rotation along parallel and perpendicular direction to the grain symmetry axis ai, r c d,ofr 
is the effective damping time due to electric dipole emission (see HDL10, HLD11), 9 is 
the angle between ai and J, and .Ftot,|| and -F tot) _i_ are total damping coefficients parallel 
and perpendicular to ai (see HDL10). In the case of fast internal relaxation, the diffusion 
coefficients A and B are averaged over the distribution function /vre 

The Langevin equation (15) is solved using the numerical integration with a constant 
timestep. At each timestep, the angular momentum Jj obtained from LEs is recorded and 
later used to find the distribution function fj with normalization / °° fjdJ = 1. 



5.2 Advantages of the Langevin Equation approach 

There are two apparent advantages of the LE approach. First, it allows us to treat the 
spinning dust emission from grains with an arbitrary grain vibrational temperature. Second, 
the impulsive excitation by single-ion collisions, which can deposit an amount of angular 
momentum greater than the grain angular momentum prior the collision, is easily included 
in Equation (14) (see [41]). Below, we briefly discuss the effect of impulsive excitations 
arising from single-ion collisions. 

DL98b showed that for grains smaller than 7 A the angular impulse due to an individual 
ion-grain collision may be comparable to the grain angular momentum prior the collision. 
Thus, infrequent hits of ions can result in the transient rotational excitation for very small 
grains. 

Let 

r icoii t> e the mean rate of ion collisions with the grain, given by 



r icoll 



/(Z g = 0)n;7ra 



J2 f( z s) n ^ a 
z g ^o 




(20) 



where $ = (2Zfe 2 /ak-BT gas ) 1 ^ 2 , g(x) = 1 — x for x < and g(x) = e x for x > 0, and /(Z g ) 
is the grain charge distribution function. The probability of the next collision occurring in 
[t, t + dt] is 



dP = TjJn exp (-t/Ti C oii) dt. 



(21) 
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The rms angular momentum per ion collision {5 J 2 ) is inferred by dividing the total rms 
angular momentum to the collision rate, and its final formula is given in Hoang et al. [36]. 

Provided that the random moment of a single-ion collision is obtained from Equation 
(21), the angular momentum that the grain acquires through each single-ion collision can 
easily be incorporated into the Langevin equation (14). Hoang et al. [36] found that the 
impulsive excitations ions extend the distribution of grain angular momentum to the region 
of high angular momentum (see next Section for its effect on spinning dust emission). 



6 Spinning Dust Emissivity 

6.1 Spinning grain of triaxial ellipsoid shape 

An irregular grain rotating with a given angular momentum J radiates at frequency modes 
u k = u) m with m = 0,±1,±2... and u k = u n with n = 1,2,3... (see Equations 12 and 13). 
For simplicity, let denote the former as u mi and the latter as u ni where i indicates the value 
for m and n. These frequency modes depend on the parameter q(s), which is determined by 
the internal thermal fluctuations within the grain. 

To find the spinning dust emissivity by a grain at an observational frequency v, first we 
need to know how much emission that is contributed by each mode oo k . 

Consider an irregular grain rotating with the angular momentum J, the probability of 
finding the emission at the angular frequency u depends on the probability of finding the 
value to such that 

pdf(u\J)du = /vrb(s, J)ds = Aexp l- J^ ) ds, (22) 

\ ^-'l^B-'vib/ 

where we assumed the VRE regime with /vre given by Equation (6). 
For the mode u = u) k {s), from Equation (22) we can derive 

pdf k {u\J) = Ms, J). (23) 

The emissivity from the mode k is calculated as 

3t, k = ^j^P ed , k (J,q<)fj(J)pdf k (u;\J)27rdJ 

~ r PcdMq>)fj(J)pdf k (u\J)2n dJ, (24) 

47T J Ji 



where q< and q > denote q < q sp and q > q sp , respectively; Ji and J u are lower and upper 
limits for J corresponding to a given angular frequency co k (J, q) = u, and 2% appears due to 
the change of variable from v to uj. 

Emissivity by a grain of size a at the observation frequency v arising from all emission 
modes is then 

£ = (25) 

k 

Consider for example the emission mode k = mo. For the case I2 slightly larger than ^3, 
this mode has the angular frequency u mo = (<fi) = (J/Ii)qo with q obtained from calculation 
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of Wnm, is independent on q for q < g S p-~ As a result 

pt^ lno HJ) = 5(w-(J// 1 )go). (26) 
Thus, the first term of Equation (24), denoted by j% mo is rewritten as 

1 r J ™ 

Ihfj(Jo) 



j; mo ,< = o /. Ped,m (J,q<)fj(J)S(uJ-(J/Ii)qo)dJ, 



2 go 



-Ped >mo (^o,g(s)), (27) 



where Jo = hu/qo, and the value of g(s) remains to be determined. 

For q > q sp , (<fi) is a function of q. Hence, the emissivity (24) for the mode k = mo 
becomes 

Ju,mo = o / « sP ed,mo(</o,g(s))jVRE(J ,s) 

2 g Jo 

+ ^ Pcd, mo ( J, q>)fj(J)pdf mo (uj\J)dJ: (28) 

where s sp is the value of s corresponding to q — q sp , and the term -P e d,m (^0; g( s )) m Equation 
(27) has been replaced by its average value over the internal thermal distribution /vre- 
The emissivity per H is obtained by integrating j° over the grain size distribution: 

3 - = -\ da^jt , (29) 
where is given by Equation (25). 

6.2 A degenerate case: grains of disk-like shape 

The spinning dust emissivity from disk-like grains (e.g., 1% = I3) is a degenerate case of 
triaxial grains. Basically, a disk-like grain with an angular momentum J radiates at four 
frequency modes: 

u mi = 4> + tip = — [h + i(l — h) cos 9] , (30) 

co ni = ip = — (1 — h)cos9, (31) 

where i = and ±1 (see HDL10; [1]). 

The emission power of these modes are given by following analytical forms (HDL10; [70]): 

2/4 A 

P ^ = 3^-0 sin0 ' ( 32 ) 
P Ua±1 = ^0,^(1 ± cos fl) a , (33) 

P. ni = 2 -§< sin29 - (34) 



2 qo approaches I1/I2 as ^3 — > I2, i.e., when irregular shape becomes spheroid. 



AME Special Issue, Thiem Hoang and A. Lazarian 



13 



For the disk-like grain, from Equation (3), the number of states in phase space s for q 
spanning from 1 — q becomes 



'h-qV 12 
h-1 



l-cos0, (35) 



where q = 1 + (h — 1) sin 2 9 has been used. Thus, for an arbitrary mode with frequency Uk, 
we obtain 

pdfk(uj\J)du = /vre(s, J)ds = /vre(0, J) sin9d9. (36) 
Taking use of u; = u)k(J,9), we derive 

pdf k (u\J) = / V re(0, J) sin#. (37) 

Therefore, by substituting Equations (32)-(34) in Equation (24), the emissivity at the 
observation frequency v = u/(2tt) from a disk-like grain of size a is now given by 

+^fV pdf mi (u\J)fAJ)dJ 



2 6c 3 Jj t 



+l^ 4 //" pdf ni (u\J)fj(J)dJ, (38) 



where pdf m±1 and pdf ni are easily derived by using Equation (37) for u m±1 and u ni , and 
J; = I\\u/(2h — 1) and J u = I\\u for m±i mode, and J\ = I\\u/(h — 1) and J u = oo for ni 
mode (see Equations 30 and 31) . 



6.3 Emissivity 

Hoang et al. [41] assumed that smallest grains of size a < 02 = 6 A, have irregular shape, 
and larger grains are spherical. To compare the emissivity from an irregular grain with that 
from a disk-like grain, they considered the simplest case of the irregular shape in which the 
circular cross-section of the disk-like grain is adjusted to the elliptical cross-section. The 
emission by two grains of different shapes with the same mass M and thickness L is under 
interest, therefore, the semi-axes of the elliptical disk is constrained by the grain mass: 

M = nR 2 L = nb 2 b 3 L, (39) 

where R = (4a 3 /3L) 1 / 2 is the radius of the disk-like grain, 62 and 63 are the length of semi- 
axes &2 and a3, and b\ = L is kept constant. Assuming that the circular disk is compressed 
by a factor a < 1 along a 2 , then Equation (39) yields 



6 2 = Q.R, 63 = a R. 



(40) 
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Figure 3: Emissivity per H from irregular grains of different degrees of irregularity 77 = 63/62 
with T vi b = 60 and 200 K in the WIM for the case in which the electric dipole moment is 
isotropically oriented in the grain body (i.e., case 1 with \i\ = [ij a/3). The emission spectrum 
shifts to higher frequency as 77 decreases (i.e., grain becomes more irregular). Here the grain 
mass is held fixed as 77 changes. Figure reproduced from Hoang et al. [41]. 



Denote the parameter rj = 63/62 = or 2 , then the degree of grain shape irregularity is 
completely characterized by 77. 

For each grain size a, the parameter 77 is increased from rj = 1 to rj = 77 max . However, 
?7max is constrained by the fact that the shortest axis &2 should not be shorter than the grain 
thickness L. The value 77 max ~ 3/2 is conservatively chosen. 

Although the irregular grain can radiate at a large number of frequency modes, only 
the modes with the order \m\ < 2 are important. The higher order modes contribute less 
than ~ 0.5% to the total emission, thus they are neglected. Hoang et al. [41] assumed that 
grains smaller than 02 have a fixed vibrational temperature T v ib (see Hoang et al. [41] for 
the detailed treatment of T vi b distribution), and that for the instantaneous value of J the 
rotational energy has a probability distribution /Vre (i-e. VRE regime, see Eq.6). 

The grain size distribution dn/da from Draine & Li [23] is adopted with the total to 
selective extinction Ry = 3.1 and the total carbon abundance per hydrogen nucleus 6c = 
5.5 x 10~ 5 in carbonaceous grains with a m i n = 3.55 A and a max = 100 A. 

The spinning dust emissivity is calculated for a so-called model A (similar to DL98b; 
HDL10), in which 25% of grains have the electric dipole moment parameter (3 = 2(3q, 50% 
have (3 = (3 and 25% have = 0.5/3o with (3 = 0.4 D. In the rest of the paper, the notation 
model A is omitted, unless stated otherwise. 

The left panel in Figure 3 shows the spinning dust emissivity for different degrees of 
irregularity 77 and with a dust temperature T v ;b = 60 K in the WIM. The emission spectrum 
for a given T v ib shifts to higher frequency as 77 decreases (i.e. the degree of grain irregularity 
increases), but their spectral profiles remain similar. The right panel shows the increase of 
peak emissivity J pea k with increasing 77. 

One particular feature in Figure 3right is that for axisymmetric grains (77 = 1), the 
emissivity increases by a factor of 1.3 with T v ib increases from 20 — 200 K. However, for the 
irregular grain with high triaxiality 77 = 1.5, the emissivity increases by a factor of 2. The 
peak frequency is increased by a factor of 1.4. 

This feature is easy to understand because the irregular grain radiate at more frequency 
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Figure 4: Emissivity per H obtained for WIM without ionic impulses using the Fokker-Planck 
equation from Ali-Haimoud et al. [1] and with impulses using our LE simulations for grain 
wobbling. The spectra are efficiently broadened as a result of impulses (see blue line). Figure 
reproduced from Hoang et al. [36]. 



modes than the axisymmetric grain. As a result, for the grain temperature increases to a 
sufficiently high value, it results in the uniform distribution of the angle between grain sym- 
metry axis and angular momentum, so that the spinning dust emissivity becomes saturated. 
On the other hand, for the triaxial grain, as T vi b increases, it allows grain to rotate about 
its axis of minimum inertia (smallest moment of inertia). As a result, the grain radiates at 
frequency modes with higher frequency, and power. 

In the case of efficient IVRET, vibrational energy is converted to rotational emission, 
which results in the increase of both emissivity and peak frequency. As shown, the energy 
transfer is more efficient for more irregular grain. The reason for that is, the more irregular 
grain allows grain spends a larger fraction of time rotating along the axis of minor inertia. 

The effect of impulsive excitations by single-ion collisions is shown in Figure 4. One can 
see that the impulses from ions can increase the emissivity by ~ 23%, and slightly increase 
the peak frequency (see Figure 4). The tail of high frequency part is obviously extended due 
to the contribution from ionic impulses with large angular momentum. 

7 Constraining spinning dust parameters and implica- 
tions 

Spinning dust emission involves a number of parameters, including grain physical parameters 
and environment parameters. Among them, the grain dipole moment and gas density are 
two most important parameters, but they can be constrained using theoretical modeling 
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Figure 5: Decrease of the peak frequency z/ pca k of spinning dust spectrum with the lower 
cutoff of grain size distribution a m i n for various environment conditions. Figure reproduced 
from Hoang et al. [41]. 



combined with observation data (see e.g., Dobler et al. [15]; Hoang et al. [41]). In the 
following, we discuss a number of parameters, which are shown to be important but more 
difficult to constrain through observation. 

7.1 Lower cutoff of grain size distribution a m i n 

The spinning dust emission spectrum is sensitive to the population of tiny dust grains, and its 
peak frequency is mostly determined by the smallest PAHs. Let a m ; n be the size of smallest 
PAHs. When a min is increased, the peak frequency t'pcak decreases accordingly. 

Figure 5 shows the variation of t'peak as a function of a m ; n for various environments, for 
the case in which the grain dipole moment lies in the grain plane (case 2) with \i\ — 
and with the VRE regime (Td = 60 K). As expected, z/ pca k decreases generically with a m ; n 
increasing. Thus, in addition to grain dipole moment, the lower cutoff of grain size also plays 
an important role. 

7.2 Constraining the shape of very small grains 

Very small grains and PAHs are expected to be nonspherical. However, contraining grain 
triaxiality using spinning dust appears rather challenging. In the simplest case where the 
grain shape can be approximated as an triaxial ellipsoid, the possibility is still low because 
there are many parameters involved in spinning dust. 



AME Special Issue, Thiem Hoang and A. Lazarian 



17 




10.0- 



i 1 1 — i — r 

model A 
case 2 



10 



100 



v (GHz) 



Figure 6: Spinning dust emissivity per H in the presence of compressible turbulence with 
sonic Mach number M s = 2 and 7, compared to that from uniform medium with nn = nn 
for the CNM (red) and WIM (green). The peak emissivity is increased, and the spectrum is 
shifted to higher frequency due to compressible turbulence. Case 2 (/ii = 0) of n orientation 
is considered. Figure reproduced from Hoang et al. [41]. 

7.3 Can compressible turbulence be observed through spinning 
dust emission? 

The discussion of interstellar conditions adopted in DL98 and other works on spinning dust 
was limited by idealized interstellar phases. It is now recognized that turbulence plays an 
important role in shaping the interstellar medium. 

For spinning dust, the turbulence can increase the emissivity due to its nonlinear de- 
pendence on material density. Indeed, in a medium with density fluctuations, the effective 
emissivity is 



where f(x)dx is the fraction of the mass with p/{p) G (x,x + dx). We use compression 
distributions f(x) obtained from MHD simulations for M s = 2 and 7 to evaluate (j u ) for the 
WIM and CNM, respectively. 

We assume the case 2 (fii = 0) of /x orientation. The resulting effective emissivity is 
compared with the emissivity from the uniform medium in Figure 6. It can be seen that 
the turbulent compression increases the emissivity, and shifts the peak to higher z/ pea k. The 
increase of emissivity is significant for strong turbulent medium. 

The distribution of phases, for instance, CNM and WNM of the ISM at high latitudes 
can be obtained from absorption lines. Similarly, studying fluctuations of emission it is pos- 




(41) 
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sible to constrain parameters of turbulence. In an idealized case of a single phase medium 
with fluctuations of density with a given characteristic size one can estimate the value of 
the 3D fluctuation by studying the 2D fluctuations of column density. A more sophisticated 
techniques of obtaining sonic Mach numbers 3 have been developed recently (see Kowal et 
al. [46]; Esquivel & Lazarian [28]; Burkhart et al. [5]). In particular, Burkhart et al. [6], 
using just column density fluctuations of the SMC, obtained a distribution of Mach numbers 
corresponding to the independent measurements obtained using Doppler shifts and absorp- 
tion data. With such an input, it is feasible to quantify the effect of turbulence in actual 
observational studies of spinning dust emission. 

7.4 Effect of dust acceleration on spinning dust emission 

Collisions of ultrasmall grains with ions and neutrals in plasma appear to be a dominant 
mechanism of rotational excitations for spinning dust emission, particularly, in dark clouds 
where UV photons are blocked out. Current spinning dust models assume Brownian motion 
of grains relative to gas, but it is known that grains may move with suprathermal velocities 
due to acceleration by turbulence (see e.g., [17]; Yan & Lazarian [77]; [78]; Hoang et al. [42]), 
and random charge fluctuations (Ivlev et al. [43]; Hoang & Lazarian [40]). The latter mech- 
anism, namely random charge fluctuations-induced-acceleration, is found to be efficient for 
tiny grains (Hoang & Lazarian [40]). 

The resonant acceleration by fast modes of MHD turbulence, which occurs when the 
grain gyroradius is comparable to the scale of turbulence eddy (i.e., r g ~ k~ l ) is considered 
a dominant mechanism for large grains (> 10~ 5 cm), whereas it is negligible for ultrasmall 
grains because the grain gyroradius falls below the cutoff scale of the turbulence due to 
viscous damping (see Yan et al. [78]; Hoang et al. [42]). 

In highly ionized media (e.g., WIM, HII regions) the resonant acceleration by MHD 
turbulence may become important for ultrasmall grains because the damping cutoff of MHD 
turbulence is suppressed due to the decrease of viscous neutral damping. We also note that 
recent observations by Paladini et al. [66]) revealed that PAHs and ultrasmall grains may be 
present in HII regions, as shown through their 8/im and 24/xm emission features, respectively. 

Thus, assuming that grain rotational kinetic energy is equal to its translational energy, 
the acceleration by these aforementioned processes is expected to increase spinning dust 
emission. Further studies should take into this issue into account. 

8 Polarization of Spinning Dust Emission and align- 
ment of ultrasmall grains 

8.1 Polarization of anomalous microwave emission 

Spinning dust emission is an important foreground component that contaminates with the 
CMB radiation in the frequency 10 — 90 GHz. The understanding of how much polarized 
this emission component is becoming an pressing question for future CMB B-mode missions. 

Recent observational studies (Dickinson et al. [14]; Lopez- Caraballo et al. [62]; Macellari 
et al. [63]) showed that the average polarization of AME is between 2 — 5%. In the last 

3 It may be seen that Alfven Mach numbers have subdominant effect on the distribution of densities (see 
Kowal et al. [46]). Thus in our study we did not vary the Alfven Mach number. 
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several years, significant progress has been made in understanding spinning dust emission, 
both theory and observation; but the principal mechanism of alignment of ultrasmall grains 
is not well understood. 

8.2 Alignment of ultrasmall dust grains 

Grain alignment is an exciting problem (see Lazarian [49] for a review). The most promising 
mechanism for the grain alignment is based on radiative torques. Proposed originally by 
Dolginov & Mitrofanov [16] it is related to the interaction of unpolarized radiation with 
irregular grains. The numerical studies in Draine & Weingartner [25] and [26] showed the 
efficiency and promise of the radiative torques (which later were termed RATs). The physical 
picture of the RAT alignment and a detailed study of important relevant effects is presented 
in Lazarian & Hoang [55], [56] and Hoang & Lazarian [37], [38], [39]. However, the efficiency 
of RATs plummets as the size of grains gets much smaller than the radiation wavelength. 
Therefore, this mechanism, which seem to provide a good correspondence with the optical 
and infrared data (see Lazarian [49], Whittet et al. [75]) cannot be applicable to ultrasmall 
spinning dust. 

Microwave emission from spinning grains is expected to be polarized if grains are aligned. 
Alignment of ultrasmall grains (essentially PAHs) is likely to be different from alignment 
of large (i.e., a > 10 -6 cm) grains discussed above. One of the mechanisms that might 
produce the alignment of the ultrasmall grains is the paramagnetic dissipation mechanism 
proposed by Davis & Greenstein [8]. The Davis- Greenstein alignment mechanism (Davis 
& Greenstein 1951, Roberge & Lazarian 1999) is straightforward: for a spinning grain the 
component of interstellar magnetic field perpendicular to the grain angular velocity varies in 
grain coordinates, resulting in time-dependent magnetization, associated energy dissipation, 
and a torque acting on the grain. As a result grains tend to rotate with angular momenta 
parallel to the interstellar magnetic field. 

Lazarian & Draine [51] (henceforth LDOO) found that the traditional picture of param- 
agnetic relaxation is incomplete, since it disregards the so-called "Barnett magnetization" 
(Landau & Lifshitz [47]). The Barnett effect-the inverse of the Einstein-De Haas effect- con- 
sists of the spontaneous magnetization of a paramagnetic body rotating in field-free space. 
This effect can be understood in terms of the lattice sharing part of its angular momen- 
tum with the spin system. Therefore the implicit assumption in Davis & Greenstein [8] 
that the magnetization within a rotating grain in a static magnetic field is equivalent to the 
magnetization within a stationary grain in a rotating magnetic field is clearly not exact. 

LDOO accounted for the "Barnett magnetization" and termed the effect of enhanced para- 
magnetic relaxation arising from grain magnetization "resonance paramagnetic relaxation" . 
It is clear from Figure 7 that resonance paramagnetic relaxation persists at the frequencies 
when the Davis-Greenstein relaxation vanishes. However the polarization is marginal for 
v > 35 GHz anyhow. The discontinuity at ~ 20 GHz is due to the assumption that smaller 
grains are planar, and larger grains are spherical. The microwave emission will be polarized 
in the plane perpendicular to magnetic field because the angular momentum is partially 
aligned with the magnetic field. 

8.3 Constraining the alignment of ultrasmall grains 

Can we constrain the alignment of ultrasmall grains through polarization of mid-infrared 
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Figure 7: Polarization for both resonance paramagnetic relaxation and Davis-Greenstein 
relaxation for grains in the cold interstellar medium as a function of frequency (from LDOO). 
For resonance relaxation the saturation effects are neglected, which means that the upper 
curves correspond to the maximal values allowed by the resonant paramagnetic mechanism. 
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(2 — 12/im) emission features? 

The answer to this question is "probably not". Indeed, as discussed earlier, mid-infrared 
emission from ultrasmall grains, takes place as they absorb UV photons. These photons 
raise grain vibrational temperature, randomizing grain axes in relation to its angular mo- 
mentum (see Lazarian & Roberge [59]). Taking values for Barnett relaxation from Lazarian 
& Draine [50], we estimate the randomization time of the 10~ 7 cm grain to be 2 x 10 -6 s, 
which is less than the grain cooling time due to IR emission. As a result, the emanating 
infrared emission will be polarized very marginally. If, however, Barnett relaxation is sup- 
pressed, the randomization time is determined by inelastic relaxation (Lazarian & Efroimsky 
[52]) and is ~ 0.1 s, which would entail a partial polarization of infrared emission. 

Can we constrain the alignment of ultrasmall grains via the ultraviolet polarization? 
PAHs and ultrasmall grains that produce spinning dust emission are likely the same particles 
that produce the prominent UV absorption feature at 2175A (see e.g., Draine & Li [24]). The 
lack of polarization excess at 2175A is consistent with the expectation that the PAHs are 
poorly aligned. However, the small degree of polarization (see Wolf et al. [76]) indicates that 
there must be some residual alignment of ultrasmall grains. The constraint for such a residual 
alignment can be obtained by fitting theoretical model with the UV polarization of starlight 
(Martin [64]). When the residual alignment available, one can predict the polarization level 
of spinning dust. 

Apart from emission from spinning dust, another new type of emission from dust is 
possible. Draine & Lazarian [21] noticed that the strongly magnetized material is capable 
of producing much more microwave thermal emission compared with non-magnetic grains. 
They suggested this as a possible alternative to spinning dust emission, which can be re- 
sponsible for a part or even most of the anomalous microwave emission. Such an emission 
can be strongly polarized, making anomalous emission an important contaminant in terms 
of CMB polarization studies. 

Further research showed that at the frequencies 20-90 GHz the spinning dust domi- 
nates. However, Draine & Hensley [18] performed new calculations of microwave response 
of strongly magnetic grains. At higher frequencies this new extensive study of evaluating 
microwave emissivity of strongly magnetic grains showed that magneto-dipole response of 
interstellar dust may be extremely important. 

9 Summary 

The principal points discussed above are as follows: 

• The model of spinning dust emission proposed by DL98 proved capable of explain- 
ing anomalous microwave emission, and its predictions were confirmed by numerous 
observations since the introduction of the model. 

• The DL98 spinning dust model has been improved recently by including the effects of 
thermal fluctuations within dust grains, impulsive excitations with single ions, transient 
heating by UV photons, triaxiality of grain shape, and compressible turbulence, which 
made the spinning dust model more realistic. 

• Spinning dust emission involves a number of grain physical parameters and environ- 
ment parameters. With the latest progress on theoretical modeling and observations, 



AME Special Issue, Thiem Hoang and A. Lazarian 



22 



the possibility of using spinning dust as a diagnostics tool for physical parameters of 
ultrasmall dust is open. 

• The spinning dust emission is expected to be partially polarized, but further studies 
on alignment of ultrasmall grains and modeling of spinning dust polarization is vitally 
required. 

AL acknowledges the support of the Center for Magnetic Self-Organization and the NASA 
grant NNX11AD32G. TH is grateful to Peter Martin for fruitful discussions on constraining 
the alignment of ultrasmall grains and polarization of anomalous microwave emission. 
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